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Abstract 

The precision measurement of the mass of the W boson is an important goal of the Fermilab 
Tevatron and the CERN Large Hadron Collider (LHC). It requires accurate theoretical calcula- 
tions which incorporate both higher-order QCD and electroweak corrections, and also provide an 
interface to parton-shower Monte Carlo programs which make it possible to realistically simulate 
experimental data. In this paper, we present a combination of the full 0(a) electroweak correc- 
tions of WGRAD2, and the next-to-leading order QCD radiative corrections to W — > Iv production 
in hadronic collisions in a single event generator based on the POWHEG framework, which is able to 
interface with the parton-shower Monte Carlo programs Pythia and Herwig. Using this new com- 
bined QCD+EW Monte Carlo program for W production, we provide numerical results for total 
cross sections and kinematic distributions of relevance to the W mass measurement at the Tevatron 
and the LHC for the processes pp,pp — > — > u^f^. In particular, we discuss the impact of EW 
corrections in the presence of QCD effects when including detector resolution effects. 



*Electronic address: C.Bernaciak@ThPhys.Uni-Heidelberg.de 
^Electronic address: dow@ubpheno.physics.buffalo.edu 



1 



I. INTRODUCTION 



The precision measurement of the mass of the H^-boson, Mw, is an important goal of the 
Fermilab Tevatron flfll and the CERN Large Hadron Collider (LHC). With precise knowl- 
edge of Mw and the top quark mass, m t , indirect information on the mass of the Higgs boson, 
Mh, within the Standard Model (SM) can be extracted from the M# dependence of radiative 
corrections to the W mass. Global SM fits to all electroweak precision data performed by 
the GFITTER collaboration and the LEP Electroweak Working Group predict the SM Higgs 



mass to be M H = 96±|i GeV 



» BE 



5| and M H = 921$ GeV (68% C.L.) |6|, respectively, 



(J ar 



which lies well within the mass range presently probed by the Tevatron [7| and LHC 
experiments. Future more precise measurements of the W and top quark masses together 
with improvements in the SM predictions of My/ are expected to considerably improve the 
indirect determination of Mh- At the Tevatron, for an integrated luminosity of C = 10 fb _1 , 



an ultimate precision of 5M W =15 MeV 
estimates range from SM W 



7 MeV 



or the W mass may be possible 101 ] . For the LHC, 



12| for C = 10 fb _1 , depend- 



ill to SM W = 20 MeV 

ing on the assumptions made for detector resolutions and theoretical uncertainties. With a 
dedicated program [13], one may be able to achieve 5Mw = 0(10 MeV). 

In hadronic collisions, the W boson mass can be determined from the transverse mass 
distribution of the lepton pair, Mt(Iv), originating from the W decay, W — > £v, and the 
transverse momentum distribution of the charged lepton or neutrino. Both QCD and elec- 
troweak (EW) corrections play an important role in the measurement of W observables at 
hadron colliders. It is imperative to control predictions for observables relevant to W pro- 
duction at least at the 1% level. Also the transverse momentum distribution of the W boson 
is an important ingredient in the current W mass measurement at the Tevatron (see, e.g. 
Ref. [10| for a review). In lowest order (LO) in perturbation theory, the W boson is produced 
without any transverse momentum. Only when QCD corrections are taken into account does 
the W boson acquire a non- negligible transverse momentum, p!f . For a detailed understand- 
ing of the distribution, it is necessary to resum the soft gluon emission terms, and to 
model non-perturbative QCD corrections. This has been done either by using calculations 
targeted specifically for resummation and parametrizing non-perturbative effects (see e.g. 
Refs. [ijl and 15]), or interfacing a calculation of W boson production at next-to-leading 
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order (NLO) in QCD with a parton-shower Monte Carlo (MC) program and tuning the pa- 
rameters used to describe the non-perturbative effects. This approach has been pursued in 
Refs. 



16 



18], for instance. Fixed higher-order predictions beyond NLO are known for fully 



differential distributions through next-to-next-to-leading order in QCD [l9M2l|. and recently 
first steps towards a calculation of the complete mixed EW-QCD 0(a s a) corrections to the 



22 



Drell-Yan production process were made in Ref. 

While QCD corrections only indirectly affect the W mass extracted from the Mt(Iv) 
distribution, EW radiative corrections can considerably distort the shape of this distribution 
in the region sensitive to the W mass. For instance, final-state photon radiation is known 
to shift M w by 0(100 MeV) [l[ [2], [23]- 28]. In the last few years, significant progress in our 



understanding of the EW corrections to W boson production in hadronic collisions has been 



made. The complete 0{a) EW radiative corrections to pp — > W ± —> t^v (£ = e, /1) were 



calculated by several groups [29H35I] and found to agree |36|, |37] . First steps towards goin. 



beyond fixed-order in QED radiative corrections in W production were taken in Refs. 



38 



421 ]. for instance, by including the effects of final-state multiple photon radiation. For a 



review o: 
Refs. 361 



the state-of-the-art of predictions for W production at hadron colliders see, e. g., 



37|, 



43[. 



As a result of all these studies, given the anticipated accuracy of a W boson mass mea- 
surement at the Tevatron and the LHC, it has become increasingly clear that it is necessary 
to not only fully understand and control the separate higher-order QCD and EW corrections, 



but also their combined effects. A first study of combined effects can be found in Ref. 
where final-state photon radiation was added to a calculation of W boson production which 
includes NLO and resummed QCD corrections. This study showed that the difference in 
the effects of EW corrections in the presence of QCD corrections and of simply adding the 
two predictions may be not negligible in view of the anticipated precision. Moreover, in 
the relevant kinematic region, i.e. around the Jacobian peak, the QCD corrections tend 
to compensate some of the effects of the EW corrections. In Ref. [45] the full set of EW 
0(a) corrections of HORACE 34J and the QCD NLO corrections to W production were com- 
bined in the MCONLO framework which is interfaced with the parton-shower MC program 



Herwig 46]. The resu 
implemented in SANC 



ts of a combination of the EW 0(a) corrections to W production as 



331 ] with Pythia [47] and Herwig can be found in Ref. [48|, without, 
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however, performing a matching of NLO QCD corrections to the parton shower. 

In this paper, we present a combination of the full EW 0(a) radiative corrections of 



Ref. 



30 



32| contained in 



the public MC code WGRAD2 and the QCD corrections to W — >• £u 



production of POWHEG-W 



17j |. One advantage of the POWHEG method 
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45 1| for the use 



in a detector simulation is that it only generates positive weighted events. Moreover, it 
provides an interface to both Herwig and Pythia. It is well suited as a starting point for 
combining EW and QCD corrections to IV-boson production in one MC program to serve 
as an analysis tool in the VT-mass measurement of the Tevatron and LHC experiments. 
The resulting MC code, called in the following POWHEG-W_EW, is publicly available at the 
POWHEG Box webpage 52( and allows the simultaneous study of the effects of both QCD 
and NLO EW corrections with both Pythia and Herwig. We do not include the effects of 
photon-induced processes and of multiple photon radiation. As has been found in earlier 



studies 



37 



42 



48], both effects, although small, still can have a non-negligible impact on 



the W^-mass measurement and should be included in view of the anticipated final precision 
of the My/ measurement at the Tevatron. This is left to a future publication. 

The technical details of our calculation and implementation of EW 0(a) corrections in 
POWHEG-W are described in Section [Til In Section II I II we first describe our crosschecks, and 
then present numerical results for total cross sections and distributions which are of interest 
for the W-mass measurement at the Tevatron and the LHC. In particular, we study the 
combined effects of EW 0(a) and QCD corrections on the Mr(/xz/ M ) and Pt(aO distributions 
in pp,pp — y W ± —> y^v^i taking into account detector resolution effects and using Pythia 
to simulate parton showering. Finally, our conclusions are presented in Section IIV1 

II. THEORETICAL FRAMEWORK 



In the following we concentrate on describing our implementation of the complete EW 
0(a) corrections to W production via the Drell-Yan mechanism q^, — y W — y //'(7) i n 
POWHEG-W. We refer to the literature for a detailed description of QCD and EW corrections to 
W boson production at hadron colliders as implemented in POWHEG-W [17J and WGRAD2 30j,|32|], 



respectively. To illustrate our implementation, we start with a schematic presentation of 



the parton-level NLO QCD cross section to W production as given in Ref. 



491 ] (see also 
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Refs. [50,1511 for a detailed description of the POWHEG BOX): 

da = £($ 2 )d$ 2 + V($ 2 )d$ 2 + [R($ 3 )d® 3 - C($ 3 )d$ 3 P] (1) 

where the 2—^3 phase space of the radiated parton is given by d& 3 = d<&2dQ va ,d, B,V,R 
denote the Born, virtual and real emission contributions, respectively, and C are counter 
terms, derived in a suitable subtraction scheme, that ensure that the term in the square 
brackets is non-singular. P denotes a projection of 2 — > 3 kinematics onto 2—7-2 kinematics. 
After some manipulation suitable for interfacing da with a parton-shower MC, the cross 
section can be written as follows |49j: 



da = B(* 2 )d* 2 [A™{0) + < LO (PT) R( l 2 ^ d) d<S> 



rad I 



(2) 



where the term in square brackets contains the Sudakov form factor A^ LO and generates the 
first emission of a light parton, while all subsequent emissions are handled by the parton- 
shower MC. B is defined as 



49] 



£($ 2 ) = B($ 2 ) + y($ 2 ) + J [i?($ 2 , $ rad ) - C($ 2 , $ ra d)]rf$rad 



(3) 



and offers a straightforward way of adding the EW corrections to the QCD 0(a s ) corrections 
contained in V and R. In detail, all changes made to POWHEG-W in order to include the 0(a) 
EW corrections of WGRAD2 are portrayed by the boxed terms contained in the _B^ 6 ($ 2 ) portion 



of POWHEG-W as follows (now we follow the notation of Ref. 
B h {$> 2 ) = [5($ 2 ) + VW$2)+|VW$2; 



5Q|): 



(4) 



a r =0 
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where the subscript 'QCD' refers to the original POWHEG-W terms. For the real terms, a r = 
corresponds to singularities occurring when the initial-state emitters are q or q' and the gluon 



could be emitted from either of them. a r — 1 corresponds to a gluon emitting an antiquark, 
q or <f , and a r = 2 to a gluon emitting a quark q' or g. The correspond to each particular 
flavor structure at the Born level where in the case of W — > iv production there are twelve. 

For each collinear piece, aw e ) = 1 corresponds to a quark/ antiquark from a hadron with 
positive (negative) rapidity emitting a collinear gluon and a©(e) = 2 to a positive(negative) 
rapidity gluon emitting a collinear quark/antiquark. 

In order to incorporate real photon emission as part of the EW 0(a) corrections, the same 
momentum used to denote the radiated parton (gluon or quark/antiquark) is used to denote 
the radiated photon. However, because our implementation of the 0(a) EW corrections 
does not include photon- induced processes, the EW contribution to the real term of Eq. @] is 
incorporated only into the a r = contribution and likewise for the collinear terms there is 
only an a eiQ = 1 term, as denoted in Eq. |U 



As described in detail in Ref. 



30| (see also Ref. |53|) we use the phase space slicing (PSS) 



method to extract the soft and collinear singular regions in the contribution of real photon 
radiation described by R^ w - In these regions the integration over the photon phase space is 
performed analytically using a soft and collinear approximation of -Re W , which is valid as long 
as the PSS parameters 5 S and 5 C are chosen to be sufficiently small. The soft part is included 
in Vew an d the remnant of the initial-state collinear singularity after mass factorization is 
denoted by G l EW , Q & y Explicit expressions for these contributions and a detailed description 
of the QED factorization scheme can be found in Ref. {30J . Finally, we refer to the appendix 
for the details of this implementation of EW 0(a) corrections into POWHEG-W. 



III. NUMERICAL RESULTS 



Numerical evaluations of the total cross sections and relevant W boson distributions have 
been obtained for both Tevatron (y/s = 1.96 TeV) and LHC (\/s = 7 TeV) processes using 
both Pythia and Herwig for QCD parton showering. After a description of the numerical 
setup in Section UlI Al and of crosschecks performed to validate the implementation of the EW 
corrections of WGRAD2 in POWHEG-W in Section IIII Bl we provide results for the total inclusive 
cross sections with and without parton showering and/or EW corrections in Section UlI CI W 
boson observables that are relevant to EW precision studies, specifically to the measurement 
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of M w , are defined and their distributions shown in Section IHIDI 



A. Setup 

The setup used to obtain the results presented in this paper closely follows Ref. [37] : 

• SM input parameters: 

• masses: M z = 91.1876 GeV, M w = 80.398 GeV, Mjj — 115 GeV 
m e = 0.51099891 MeV, = 0.1056583668 GeV, m T = 1.77684 GeV 
m u = 0.06983 keV, m c = 1.2 GeV, m t = 171.2 GeV 

m d = 0.06984 keV, m s = 0.15 GeV, m b = 4.6 GeV 

• W width: Y w = 2.141 GeV 

• EW coupling parameters: a(0) = 1/137.035999679 
cos^ = Mw/Mz, sin 2 9 W = 1 — cos 2 9 W 

• CKM matrix elements: \V ud \ = \ V CS \ = 0.975, \V US \ = \V cd \ = 0.222, 

1^1 = 1^1 = 1^1 = 1^.1 = 0,1^1 = 1 

• WGRAD2 flags: qnonr=0, QED=4, lfc=l 

• renormalization/factorization scales: \ip = /x^ = \x = W boson invariant mass 

(/^ = ^QED = ^QCD) 



Parton Distribution Function (PDF): CTEQ10 [54 1 



• bare acceptance cuts: pr{l) > 25 GeV, Pt{vi) > 25 GeV, and \r)e\ < 1 

• calometric setup: in addition to bare acceptance cuts, smearing of the four-momenta 
is applied, and we limit the photon energy for small muon-photon angles as described 
below. 

• Pythia settings: MSTP(61)=1, MSTP(71) = 1, MSTJ(41)=1 which corresponds to 
all QED showering turned off 

The calometric setup includes smearing of the final-state four-momenta to take into ac- 
count the uncertainty in the energy measurement in the detector. Gaussian smearing of 



the final-state four-momenta is simulated with DO or ATLAS inspired smearing routines. 
All observables are then calculated from the smeared momenta. Muons are detected in the 
muon chamber and the requirement that the associated track is consistent with a minimum 
ionizing particle. Therefore, for muons at the Tevatron and the LHC, we require a small 
photon energy for small muon-photon opening angles, i.e. we require that < 2 GeV for 
AR^ < 0.1 and E 1 < O.IE^ for 0.1 < AJ? M7 < 0.4. AR^ denotes the separation of a 
charged lepton and photon in the pseudo-rapidity azimuthal angle plane defined as: 



AR h = yj A<f)f 7 + A v f 7 (5) 

The results in this paper are obtained in the constant-width scheme and by using the fine 
structure constant, a(0), in both the LO and NLO EW calculation of the W observables. 
Since QED radiation has the dominant effect on observables relevant to the W mass mea- 
surement, we only include resonant weak corrections (qnonr=0), i.e. we neglect weak box 
diagrams. Their impact is important in kinematic distributions away from the resonance re- 
gion and can be studied by choosing qnonr=l. We include the full set of QED contributions 
(QED=4), i.e. initial-state and final-state radiation as well as interference contributions. 
The QED and QCD factorization and QCD renormalization scales are chosen to be equal 
and we assume that the factorization of the photonic initial-state quark mass singularities is 
done in the QED DIS scheme (lfc=l). The QED MS scheme is implemented as well (lfc=0) 
and both schemes are defined in analogy to the corresponding QCD factorization schemes. 
A description of the QED factorization scheme as implemented in P0WHEG-W_EW can be found 



in Ref. [30]. 



The fermion masses only contribute to the EW gauge boson self-energies and as regulators 
of the collinear singularity. The mass of the charged lepton is included in the phase space 
generation of the final-state four-momenta and serves as a regulator of the singularity arising 
from collinear photon radiation off the charged lepton. Thus, no collinear cut needs to be 
applied (collcut=0 in P0WHEG-W_EW) on final-state photon radiation, allowing the study of 
finite lepton-mass effects. Note that the application of a collinear cut on final-state photon 
radiation (collcut=l) is only allowed in the electron case and only when a recombination of 
the electron and photon momenta is performed in the collinear region (usually defined by 



30] for a detailed discussion). In this paper we present results for the 



Ai? e7 < Rem, see Ref. 
pp,pp —> W ± —> n^Vu, processes in both the bare and calometric setup. 



B. Crosschecks 



In order to be sure that the EW corrections are properly implemented, a number of 
crosschecks were performed. In the first, the QCD corrections in POWHEG-W_EW were turned 
off (i.e. all terms labeled with the subscript 'QCD' in Eq. H] were set to zero) and the 
numerical results of each piece of the NLO EW corrections (i.e. Vew, ^ ew / e -Rew of Eq. @| 
were compared to WGRAD2. We also compared results for the total inclusive cross section 
and the Mt(W) and px{£) distributions obtained with POWHEG-W_EW when only including 
EW 0(a) corrections with those obtained with WGRAD2. In all these comparisons we found 
good agreement within the statistical uncertainties of the numerical integration (see also 
Sections IIIICfHIDp . This is the primary indication that the NLO EW corrections were 
implemented properly using the numerical phase space integration of the POWHEG BOX. In 
the second type of crosscheck, the numerical cancellation of the PSS parameters 8 S and S c 
was tested by running POWHEG-W_EW without the NLO QCD corrections (or parton showering 
capabilities) of POWHEG-W for different choices of these parameters and observing that the 
cross sections agree within the statistical error of the numerical integration as long as the 
PSS parameters are chosen small enough so that the soft/collinear approximation is valid. 
To illustrate this cancellation we show in Table [I] the results for the total inclusive cross 
sections for W ± — >■ p^v^ production at the Tevatron and the LHC. Finally, we also checked 
that the QCD NLO cross sections still coincide with those obtained with the original code, 
POWHEG-W. 
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levatron 


LHC 


{ Si si \ 




W 


W 


0.01,0.005 


362.4(2) 


1059.0(5) 


758.7(8) 


0.01,0.001 


362.4(2) 


1059.1(7) 


759.2(5) 


0.001,0.0005 


362.3(2) 


1059.4(9) 


759.4(5) 


0.001,0.0001 


362.3(2) 


1059.2(8) 


759.3(5) 



TABLE I: Results for the total cross sections (in pb) to pp,pp — > W — > fu , for different choices 
of S s and S c parameters at the Tevatron (\/5* = 1.96 TeV) and LHC (^/S = 7 TeV). These results 
reflect the exact (weighted) NLO results of P0WHEG-W_EW when only including EW corrections, with 
bare cuts. 

C. Total Cross Sections 

In Tables [III and III II we present results obtained with P0WHEG-W_EW for the total cross 
sections of pp,pp — > W — > p^v^ processes at the Tevatron with center-of-mass (CM) 
energy v^S* = 1.96 TeV and the LHC with v^? = 7 TeV for the bare and calometric setup, 
respectively. We show results separately for the NLO EW and NLO QCD cross sections and 
combined NLO EW+QCD results with and without QCD parton shower for both Pythia 
and Herwig. The NLO EW and QCD results coincide with those that can be obtained with 
WGRAD2 and P0WHEG-W, respectively. At the level of the total cross sections, the combined 
results can be approximated by simply adding the QCD and EW cross sections. As we will 
see in Section IIIIDl this is not necessarily the case when studying kinematic distributions 
after applying parton showering. For instance, for W + production at the Tevatron(LHC) 
with bare cuts (see Table [TT]) the EW 0(a) corrections increase the LO total cross section by 
3.6%(3.4%), the combined (QCD+EW) corrections increase the QCD cross section at NLO 
by 3.3%(3.5%) and when parton-showering with Pythia is included by 3.7 ± 0.4%(3.7 ± 
0.4%). When considering the calometric setup (see Table IHip the impact of the NLO EW 
corrections is considerably reduced, and the EW 0(a) corrections increase the LO total 
cross section by only 1.6%(1.4%), the combined (QCD+EW) corrections increase the QCD 
cross section at NLO by 1.1%(1.1%) and when parton-showering with Pythia is included by 
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1.6 ± 0.3%(0.9 ± 0.3%). 





Tevatron 


LHC 


W + 


W+ 


W~ 


LO 


349.81(2) [349.77(1)] 


1024.0(1)[1023.9(1)] 


731.69(6) [731.63(2)] 


NLO EW 


362.4(2)[362.55(2)] 


1059.0(5)[1059.6(1)] 


758.7(8)[759.26(3)] 


NLO QCD 


384.66(4) 


1022.7(2) 


750.8(1) 


NLO (QCD+EW) 


397.2(2) 


1058.0(6) 


778(1) 




Pythia 


Herwig 


Pythia 


Herwig 


Pythia 


Herwig 


LO0PS 


308.4(7) 


311.8(7) 


854(3) 


866(3) 


634(2) 


639(2) 


NLO QCDigiPS 


375.3(8) 


378.5(8) 


1014(3) 


1027(3) 


744(2) 


750(2) 


NLO (QCD+EW)®PS 


389.3(8) 


393.1(8) 


1052(3) 


1066(3) 


766(2) 


774(2) 



TABLE II: Total cross section results (in pb) of P0WHEG-W_EW for — > [i^v^ production at the 
Tevatron (VS = 1.96 TeV) and the LHC (VS = 7 TeV) with bare acceptance criteria as listed 
in Section |fflAl Shown are results for the LO, NLO EW, NLO QCD, the combined NLO QCD 
and EW cross sections, as well as results including showering performed with Pythia or Herwig as 
provided by the P0WHEG Box. The errors shown in parenthesis are statistical errors of the Monte 
Carlo integration. As a cross check the results of WGRAD2 are provided as well (in square brackets). 
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levatron 


LHC 




TT/+ 

w ^ 


TT/ — 
W 


LU 


320.66(2) 


985.62(9) 


710.52(6) 


att f \ riTTr 

JNLU bw 


325.9(1) 


999.0(5) 


715.3(8) 


ATT z' \ f \f U \ 

JNLU OUD 


369.75(4) 


1037.9(2) 


758.1(1) 


ATT ( \ [ i \f * \ \ I TT'ASA 

JNLU (C4CU+tvWJ 


373.9(1) 


1049.0(6) 


765(1) 




Pythia 


Herwig 


Pythia 


Herwig 


Pythia 


Herwig 


LO^PS 


296.3(6) 


298.2(6) 


853(3) 


861(3) 


622(2) 


626(2) 


NLO QCD^PS 


358.4(8) 


361.4(8) 


1006(3) 


1019(3) 


736(2) 


743(2) 


NLO (QCD+EW)®PS 


364.1(8) 


366.9(8) 


1015(3) 


1026(3) 


743(2) 


751(2) 



TABLE III: Total cross section results (in pb) of P0WHEG-W_EW for W ± -> // ± f M production at the 
Tevatron (y/S = 1.96 TeV) and the LHC (y/S = 7 TeV) with calometric acceptance criteria as 
listed in Section MEM Shown are results for LO, NLO EW, NLO QCD, the combined NLO QCD 
and EW cross sections, as well as results including showering performed with Pythia or Herwig as 
provided by the P0WHEG Box. The errors shown in parenthesis are statistical errors of the Monte 
Carlo integration. 



D. Transverse W mass and charged lepton momentum distributions 

Differential distributions for the following VF-boson observables are shown: transverse 
mass of the W, Mt(W), and the transverse momentum of the muon, Pt{h)- The transverse 
W mass is defined in terms of lepton-neutrino pair observables as 

M T (W) = V2M^tH(1-cos(A<M,)) (6) 



with <j) the azimuthal angle of the charged lepton or neutrino and the difference between 
them. Both observables are being used to perform a high precision W mass measurement 
at the Tevatron jl, 2|. The W mass extracted from these observables is especially sensitive 
to changes in the lineshape in the vicinity of the Jacobian peak. As has been well-studied 
in the literature, final-state photon radiation greatly affects the distributions in this region 
and predictions for these effects need to be under good theoretical control. Here we will not 
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provide a detailed phenomenological study of these EW effects, which are available in the 
literature (see, e.g., Ref. {37I for an overview), but rather shall explore how the impact of 
EW 0(a) corrections is affected by the presence of QCD radiation when considering realistic 
lepton identification criteria. We only briefly illustrate the main features of the impact of 
the EW 0(a) corrections on the Mt(W) distribution in Fig. [1] and on the pt{1 = yU, e) 
distribution in Fig. |5]when considering the bare and calometric setup by showing the relative 
corrections defined as 

W%) = dC \^ d ° x 100 . (7) 

dO 

The large distortion of the Jacobian peak, especially when only bare cuts are applied, is due 
to collinear final-state photon radiation which results into large logarithmic enhancements 
of the form a log(m|/s) , where mi denotes the charged lepton mass and s the partonic 
CM energy squared. In the electron case, when more realistic experimental conditions are 
taken into account, the electron and photon four-momentum vectors are recombined to an 
effective electron four-momentum vector if their separation AR ei in the azimuthal angle- 
pseudorapidity plane is smaller than a critical value R cu t = 0.1. In that case, these mass- 
singular logarithms cancel, and only a small effect of the EW corrections survives. In the 
muon case, however, the muon is well separated from the photon so that in the calometric 
setup, the distortion of the Mt{W) and Pt{h) distributions around the Jacobian peak is 
more pronounced than in the electron case. Since the electron case is very sensitive to the 
details of the lepton identification requirements, a detailed study of combined EW and QCD 
effects should be performed in collaboration with experimentalists involved in the W mass 
measurement. We leave such a study to a later publication and concentrate in this paper on 
discussing the impact of QCD corrections on the EW effects in the muon case. Note that the 
results presented in Ref. [44j are obtained for the electron case with bare cuts, and therefore 
larger effects have been observed than it can be the case with a more realistic treatment. 
Finally, we will only present results for the W + — > process, since even at the LHC the 
EW effects in the Mt(W) and Pt{h) distributions are similar in W + and W~ production, at 
least in the kinematic region of interest (note that there is no distinction between W + and 
W~ production at the Tevatron because of the symmetric initial state). Also in the presence 
of QCD radiation we found that the relative corrections discussed in this paper exhibit very 
similar features for W + and W~ observables at the LHC. 
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FIG. 1: Relative corrections 5ew to the Mt(W) distributions for pp — > W + e + z/ e , lii^v^ at 
yfS =1.96 TeV on the left and pp — >■ W + — > e + ^ e , at yS = 7 TeV on the right obtained with 

WGRAD2. Bare and calometric cuts are shown for comparison. 



muon,bare 
muon,calo 
electron, bare 
electron,calo 





Pt(I) [GeV] 



MO [^V] 



FIG. 2: Relative corrections 5ew to pr(0 distributions for pp — > W + — > e + v e , at y/S = 1.96 
TeV, on the left and pp -> W + -> e + v e ,^ + u^ at \/S = 7 TeV on the right obtained with WGRAD2. 
Bare and calometric cuts are shown for comparison. 



In Figs. El HI E and Figs. El HOI HH [U we provide respectively M T (W) and p T (fi) 
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distributions calculated at LO®PS, NLO EW, NLO QCD®PS and NLO (QCD+EW)®PS 
as provided by POWHEG-W_EW, where we used Pythia to shower the events. We also produced 
these distributions with Herwig, but since the effects we are interested in, i. e. the change 
of the relative impact of EW corrections in the presence of QCD radiation are similar we 
only show results obtained with Pythia. We note that a detailed, tuned comparison of the 
impact of these two parton-shower MCs on W observables in the presence of the complete 
EW 0(a) corrections, which is especially important in determining their contribution to the 
theoretical uncertainty in the extraction of the W mass, can now be conveniently performed 
with POWHEG-W_EW. To illustrate the impact of the various higher-order corrections on Mt(W) 
and Pt(h) distributions we also show various relative corrections. Apart from <5 EW of Eq. El 
we show the impact of QCD corrections (<5qcd) and combined QCD and EW corrections 
(<5qcdew) hi the presence of a QCD parton shower relative to the LO parton shower result, 
defined as 



^°"QCD®PS rfoLO®PS 



<W%) = -^h^r^ >< 100 , (8) 



dO 

and 



^°"(QCD+EW)8iPa _ tfo- L Q, p> 



Ww(%) = d ° dm d ° x 100 , (9) 

dO 

Note that the EW NLO results shown in Figs. El SI El Eland Figs. El EH EH 03 can be directly 
compared to and should agree with the results obtained with WGRAD2 shown in Figs. [Hand El 
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m T (W), GeV mT (w), GeV 



FIG. 3: Mt(W) distributions and relative corrections <5ew> <5qcD) ^qcdew for pp — > W + — > 

V~S = 1.96 TeV, obtained with POWHEG-W_EW, with bare cuts. Parton showering (denoted by PS) is 

performed by interfacing with Pythia. 




m T (W), GeV mT ( W)i QeV 



FIG. 4: Mt(W) distributions and relative corrections <5ew> <5qcD) ^qcdew for pp — > W + — > 

\/~S = 7 TeV, obtained with POWHEG-W_EW, with bare cuts. Parton showering (denoted by PS) is 

performed by interfacing with Pythia. 
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m T (W), GeV m T (\N), GeV 



FIG. 5: Mt(W) distributions and relative corrections <5ew> <5qcD) ^qcdew for pp — > W + — > 

V~S = 1.96 TeV, obtained with POWHEG-W_EW, with calometric cuts. Parton showering (denoted by 

PS) is performed by interfacing with Pythia. 




m T (W), GeV m T (\N), GeV 



FIG. 6: Mt(W) distributions and relative corrections <5eWj <5qcd> ^qcdew fo r PP W + — > 

\/~S = 7 TeV, obtained with POWHEG-W_EW, with calometric cuts. Parton showering (denoted by PS) 

is performed by interfacing with Pythia. 



Overall, we observe that the impact on the Mt{W) distribution of the EW corrections, 
<5ew> alone is seen to be slightly negative in the peak region, while the effect of the QCD (NLO 
and PS), 5qcd, corrections alone consistently increases the LO®PS Mt(W) distribution, as 
expected. Now, the combined relative NLO (QCD+EW)(g>PS corrections, £qcdew, still 
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follow in magnitude the effect one expects when simply adding NLO QCDcgRS and NLO 
EW corrections, but seems to be slightly different in shape in the peak region at the LHC 
when applying bare cuts. The dip observed in the relative #ew correction seems now to be 
somewhat washed out. 

To study this possible effect of combining EW and QCD corrections more closely, we 
follow the discussion of Ref. j^j], and for each of the observables, we define 

cfcrNLO EW 

= ( 10 ) 

dO 

rfCT (QCD+EW)0PS 

r 2 = — 

dO 

and show in Figs. [S] r 12 together with their ratio 1Z = r 2 /ri. We see that within the 
statistical uncertainty of the numerical integration, TZ is largely consistent with unity, i. e. 
generally the NLO EW corrections in the presence of QCD effects behave like those of EW 
alone, only at the LHC there seems to be a slight change in shape above the peak region, but 
given the large fluctuation in this region this is more likely a relic of the numerical integration. 
Since QCD radiation is known not to have a significant effect on the shape of the Mt{W) 
distribution, it is no surprise that the main features of the EW corrections in the presence 



of QCD corrections are largely unchanged. This has also been observed in Refs. [44j, |45|. 
Our results cannot be directly compared with earlier studies such as those presented in 
Refs. [44 



45j , since in Ref. [44j results are provided for the electron case with bare cuts and 
in Ref. [45] slightly different selection criteria have been used and the EW input scheme 
has been adopted 1 . Nevertheless, the overall features of the effects of combining QCD and 
EW corrections turn out to be similar. It will be interesting to perform a tuned comparison 
of different implementations of EW/QED corrections in NLO QCD+resummed calculations 



of W boson observables, i. e. of POWHEG-W_EW, MC@NL0/H0RACE 
is left to a future publication. 



45j and ResBos-A [44j], which 



1 The results obtained in the G M scheme can be estimated from our results in the a(0) scheme by 
= (<j% LO - 2Ar<jL°)a 2 G Ja{0) 2 with Ar = 0.0315, a Gft = V^G„ sin 2 e w M^/ir, and G„ = 
1.16639- 10- 5 GcV~ 2 . 
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FIG. 7: Ratios and their ratio 1Z for pp — > W + — > \/~S = 1.96 TeV, obtained with 

POWHEG-W_EW, with bare (left) and calometric (right) cuts. Parton showering is performed by inter- 
facing with Pythia. 




FIG. 8: Ratios r\£ and their ratio TZ for pp W + — > ^ + v^-, = 7 TeV, obtained with 
POWHEG-W_EW, with bare (left) and calometric (right) cuts. Parton showering is performed by inter- 
facing with Pythia. 



We now turn to the discussion of the Pt{h) distributions, and show in the right-hand 
plots in Figs. l9lfT2l again the relative corrections <5ew of Eq. [7J <5qcd °f Eq. M and 5qcdew 
of Eq. [91 In contrast to the Mt(W) distribution, now the effects of QCD radiation on the 
shape of the Pt(aO distributions are quite pronounced, as can be seen by comparing the 
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exact NLO EW distribution with the distributions obtained with showered events shown on 
the left-hand side of Figs. I9HT21 So, we expect to see a change in how the EW corrections 
impact the Pt{^) distributions due to a non-trivial interplay of QCD and EW corrections 
in the combined result. After all, as can be seen in Eq. [2] with B of Eq. HI higher-order 
QCD-EW interference terms are present and may have a non-negligible impact on the Pt(aO 
distributions. Therefore, it is interesting to note the overall significant shape change between 
the exact NLO EW and the combined NLO (QCD+EW)<g>PS p T (fx) distributions. That the 
NLO EW corrections, in the presence of QCD effects, tend to be quite subdued is obvious in 
the magnitude and shape of the relative NLO (QCD+EW)(g>PS corrections compared to the 
NLO EW relative corrections, as shown on the right-hand side of these figures. Especially 
the dip in the Jacobian region from the NLO EW corrections with bare cuts is now mostly 
washed out. 




Pt('). GeV p T (l), GeV 



FIG. 9: Pt{^) distributions and relative corrections <5ew> <5qcd> <5qcdew f° r PP ~^ W + — > A t+ f /1 , 
V~S = 1.96 TeV, obtained with POWHEG-W_EW, with bare cuts. Parton showering (denoted by PS) is 
performed by interfacing with Pythia. 
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FIG. 10: PrifJ-) distributions and relative corrections 5ew> <5qcd> <5qcdew for pp — > W + — > 

\/~S = 7 TeV, obtained with P0WHEG-W_EW, with bare cuts. Parton showering (denoted by PS) is 

performed by interfacing with Pythia. 
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FIG. 11: Pr(fj) distributions and relative corrections <5ew> $qcd, <5qcdew for PP — > W + — > 

V~S = 1.96 TeV, obtained with P0WHEG-W_EW, with calometric cuts. Parton showering (denoted by 

PS) is performed by interfacing with Pythia. 
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FIG. 12: Pr(fJ-) distributions and relative corrections Sew, <5qcd> <5qcdew f° r PP ~^ W + — > 

y/S = 7 TeV, obtained with P0WHEG-W_EW, with calometric cuts. Parton showering (denoted by PS) 

is performed by interfacing with Pythia. 



Again, we study these effects more closely in Figs. [13] and [TH by comparing the impact of 
NLO EW corrections alone and their impact in the presence of QCD corrections as described 
by r% and r 2 of Eq. [10] respectively. Clearly, unlike in the case of the Mt(W) distribution, r x 
and T2 are now quite different so that their ratio TZ exhibits an interesting shape especially 
around the Jacobian peak. The same feature has also been observed in the 1Z ratio of 
Ref. [441 ] . As expected, in the calometric setup, this difference is much less pronounced due 
to the smaller impact of the EW corrections on the shape of the Pt{h) distributions. 
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FIG. 13: Ratios ri 2 and their ratio 1Z for pp — > W + — > = 1.96 TeV, obtained with 

POWHEG-W_EW, with bare (left) and calometric (right) cuts. Parton showering is performed by inter- 
facing with Pythia. 
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FIG. 14: Ratios ri 2 and their ratio 7Z for pp — > W + — > fx + fuj V~S = 7 TeV, obtained with 
P0WHEG-W_EW, with bare (left) and calometric (right) cuts. Parton showering is performed by inter- 
facing with Pythia. 



IV. CONCLUSIONS 

In this paper we described the combination of the complete EW O{o) corrections with 
NLO QCD+resummed corrections to W production in hadronic collisions, based on im- 
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plementing the EW corrections of WGRAD2 in POWHEG-W. Using the resulting MC program, 
POWHEG-W_EW, which is publicly available on the webpage of the POWHEG BOX, we presented 
results for the transverse W mass and charged lepton momentum distributions, taking into 
account lepton identification requirements which are closely modeled after those used in the 
high-precision measurement of the W mass at the Tevatron. In view of the anticipated 
precision of the W mass measurement at the Tevatron and the LHC, predictions for these 
observables including higher-order radiative corrections have to be under excellent control. 
Tools such as POWHEG-W_EW that allow the study of combined EW and QCD corrections are 
important in reducing the theoretical uncertainty in the W mass measurement. We espe- 
cially concentrated on studying whether there is a change of the impact of EW corrections 
when QCD radiation is present as described by POWHEG+Pythia in the kinematic region 
where EW corrections are known to have a significant impact on the extracted W mass. 
We found interesting QCD-EW interference effects in the Pt(/-0 distributions, i. e. effects 
that go beyond simply adding QCD and NLO EW corrections, that change the shape of 
the distribution around the Jacobian peak. These effects are similar to those observed in 
Refs. 44J, |45[, and their impact on the W mass extracted from the pr(/-0 distribution should 
be studied in more detail by using realistic detector resolution effects, ideally in close col- 
laboration with the experimentalists performing the W^-mass measurement. Moreover, these 
findings also suggest that a calculation of the complete mixed EW-QCD 0(aa s ) corrections 
is desirable to further reduce the theoretical uncertainty and to obtain an accurate estimate 
of the theory uncertainty due to missing higher-order corrections. Further improvements 
that are planned for POWHEG-W_EW include the implementation of the known higher-order 
QED and EW effects, i. e. beyond NLO, of photon-induced processes, and the usage of 
an updated PDF that fully considers QED corrections, once available. Since POWHEG-W_EW 
interfaces to both Pythia and Herwig it is also a convenient tool to perform a tuned compar- 
ison of QCD+EW effects when using either parton-shower MC. Finally, since POWHEG-W_EW 
includes the complete NLO EW corrections, it is interesting to note that the effects of EW 
Sudakov logarithms that become numerically important in distributions at high energies can 
now also be studied in the presence of QCD radiation. This is especially interesting for the 
search for W bosons at the LHC. 
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Note Added 

Shortly after the submission of our arXiv paper, another implementation of EW correc- 
tions to single W production into the POWHEG Box became available 55]. 
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V. APPENDIX 



In the following we describe in detail the implementation of each EW piece framed in 
Eq. HI which are added to the QCD contribution after testing if the user requests to include 
EW corrections (wgrad2=l). An attempt is made to describe both the analytical pieces in 
the code as well as their variable names. This is important so the user can be aware of how 
certain flags affect their values as well as of which part of the EW corrections is contained in 



each piece. For more explicit expressions of the EW 0(a) corrections see Refs. [29l. l30l. 132 ] . 



A. The virtual+soft finite piece, Vj^v^) 

The quantity V^($ 2 ) is defmed in POWHEG-BOX/W_EW-BW/virtual_EW. f . There are two 
flags which affect the outcome, both can be set in the proper POWHEG input file. Except for 
qnonr, all of the flags are used in the definitions of the other boxed terms in Eq. HJ For 
completeness, the flags and their descriptions are listed below. 
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• QED: This flag toggles between different subsets of NLO QED contributions. 

• QED=1 Initial State (IS) photon radiation. 

• QED=2 Final State (FS) photon radiation. 

• QED=3 interference between IS and FS photon radiation. 

• QED=4 IS, FS and interference 

Note that the gauge invariant separation into IS and FS QED contributions has been per- 



formed according to Ref. [29 ]. 



• qnonr: This flag toggles between the treatment of EW corrections to resonant/non- 
resonant W production. In the case of resonant W production (qnonr=0), a gauge 
invariant separation into weak and QED 0(a) corrections and into IS and FS con- 
tributions is available according to Ref. [29[. In this case, no weak box diagrams are 
included and the weak form factors are evaluated at s — M^r. 

• qnonr=0 excludes the weak box diagrams, corresponds to resonant W production 
only. 

• qnonr=l includes the weak box diagrams, i.e. includes the complete EW G(a) 
corrections to W production. 

In choosing qnonr=l, the full set of virtual diagrams are used, including the weak box 
diagrams which are necessary to describe non-resonant W production. Since in case of non- 
resonant W production the gauge invariant separation in IS and FS contributions according 
to Ref. 2^1 is no longer possible, the entire set of real radiation diagrams must be included. 
Therefore, if choosing qnonr=l the user must also set QED=4. 

The contents of ^^($2), without the PDF factors, are now explicitly shown. 

^($2) ~ soft (1:12) + virt(l:12) (11) 

The choice of qnonr affects virt (1 : 12) as follows. For qnonr=0 

virt (1:12) =br_born(l:12)fce[fvwp(l) IS +fvwp(2) FS ] — (12) 
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br_born(l : 12) are POWHEG-W's definition of the Born contribution and are equivalent to 
^~^2\Mb\ 2 for each flavor structure, f vwp (1 : 2) are the modified weak one-loop contributions 
for resonant W production and are defined in /POWHEG-BOX/W_EW-BW/libweak.f . The user 
has the choice to include separately the IS and FS contributions through the proper choice 
of QED. Shown above is the outcome corresponding to QED=4. For qnonr=l, 

7r 2 a 2 1 1 1 



virt(l:6) = CKM(1:6) 



|2 



4, \s- M 2 v + iT w M w \ 2 2s 3 

xa2qqw(s, -2k q ■ k h -2k q ■ k h 16(fc ? • k T ) 2 ) (13) 



virt(7:12) = CKM(7 : 12) 2 — =— = 

1 1 s A w \s- M 2 v + iT w M w \ 2 2s 3 



xa2qqw(s, -2k q ■ k h -2k q ■ k h 15{k q ■ kj) 2 ) (14) 

The function a2qqw(. . .) corresponds to the full weak one-loop contribution and is also 
defined in /POWHEG-BOX/W_EW-BW/libweak. f . The soft contributions of Eq. [TT] are given 
by 



soft(l:6) = |CKM(1:6) 



2 n 2 a 2 1 16 (kg ■ kj 



s%? 2s 3 \s-M^ + iT w M w 



2 



x [k_fac(2)i S + k_fac(2) FS + k_fac(2) int ] (15) 



soft (7: 12) = |CKM(7:12) 



|2 n 2 a 2 1 16 (k q ■ kjf 



4 2s 3 \s - M 2 , + iY w M w \ 2 



x [k_fac(l)i S + k_fac(l) FS + k_fac(l) int ] (16) 

The k_f ac(l :2) terms correspond to the finite soft photon contributions and are defined in 
/POWHEG-BOX/W_EW-BW/libqed.f . Again, the case of QED=4 is shown above, but the user 
can include any or all of these k_f ac terms by adjusting this flag. The k_f ac, and hence the 
soft contribution depends on the soft PSS parameter, S s , and when collcut=l is chosen also 
on S c . 
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1 f 1 f 

B. The collinear remnants, G E '^ , Gg'w e 

These are defined in /POWHEG-BOX/W_EW-BW/collinear_EW.f . There is also the inclusion 
of a new flag, If c: 

• If c: This flag controls the choice of the QED factorization scheme. 

• If c=0 MS scheme 

• If c=l DIS scheme 

POWHEG only implements the QCD MS scheme. Hence, as this flag only applies to the EW 

portions, one has the option to choose different schemes for the QED and QCD factorization. 

We now describe analytically the contents of these contributions. 

_ii.fi ^vi i-6 , , x^ 71 " 2 ^ 2 1 16 „ , N o kn_iacborn 

+ G^e = ICKMd :6) f-j- -- (k, . k f ^L——^ 



1_ f x 1 \ splitzl \ i — / v f%2 \ splitz2 

-qi — ,Hr \ q 2 {x 2 ,ii R ) (1 - xx) + -qi{x 1 ,fx R )q 2 — (1 - x 2 ) 

9 Ui / zi 9 \z 2 J z 2 



-v ' s - 



^©,EW u e,EW 



^vi 7-i2 ^,17-12 , , x,9 7r2 « 2 1 16 , ,o kn_iacborn 



(17) 



X 



4 /xi \ splitzl 1 /ac 2 \ splitz2 

-qi — q 2 {x 2 ,ii R ) (1 - Xi) + -q 1 {x 1 ,fx R )q 2 — (1 - x 2 j 

9 Vzi / Zi 9 \z 2 / z 2 



, r 1 ^r 1 

U ®,BW u 0,EW 



splitzl and splitz2 are functions of z% and z 2 respectively. For a generic z they are 



a 

splitz = — 

2rr 



1 + z , / s z 6 C . 

In -=- 7 — — + 1 - ^ - If c * f collz 



i-z V/4(i--2) 22 



with 

1 + z 2 , (l-z\ 3 



(19) 



f collz = _ In ^— j - ^) + + 3 ( 2 °) 
In Eq. [19] above, one sees explicitly the QED factorization scheme dependence. Note that 
the soft part of the QED PDF counterterm is included in k_f ac. Also, the dependence on 
the PSS parameter S c is explicitly shown. 
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C. The real contribution, R 



EW 



R^ w are defined in /POWHEG-BOX/W_EW-BW/real_EW.f . As mentioned earlier, POWHEG-W 
considers not only qq, but gq and gq induced processes, while the calculation of WGRAD2 only 
includes qq induced processes. This is why in Eq. H] there is only one i?g W contribution 
corresponding to a r = 0. Also, because the EW calculation was performed with two-cutoff 
phase-space-slicing this EW contribution is finite. Therefore, care is taken to only return 
values of R^ w which pass certain criteria, namely that they be away from the soft or collinear 
regions of the photon phase space. One can also consider subsets of the real corrections by 
adjusting the flag QED when qnonr=0. The real EW contributions, up to the PDF factors, 
are proportional to subsets of the QED radiation matrix element squared as follows 



fjj Sp,C0l : ' ' 



fb sp,col 

EE 

f b sp,col 



\Mjs\1 + \M FS \ 2 fb + 2Ke [M ls M* s ] h 

QED=1 QED=2 



QED=3 



(21) 



(22) 



QED=4 



The integration of R^w over ^ ne photon phase space is finite after applying soft and collinear 
cuts as shown in Eq. HI The dependence on the PSS parameters is canceled numerically be- 
tween the contributions describing soft, collinear and real hard photon radiation as discussed 
in Section 1111 Bl The default values set in /POWHEG-BOX/W_EW-BW/virtual_EW . f represent 
an optimal choice and should only be changed with care. 
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